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1. Introduction 

In a recent article pQ, we presented variational calculations of phase shifts for the 
elastic scattering of very low energy positrons, e + , by ground state molecular hydrogen, 
H2. Those calculations involved modifications of the variational method originally due 
to Kohn [2] and Hulthen [31 H] (see also Rubinow [5j), which is the analogue of the 
Rayleigh-Ritz method [6] widely used in variational calculations on bound states. Our 
earlier analysis concentrated on the so-called Schwartz singularities [8] , which are well 
known to be associated with anomalous behaviour in the results of Kohn calculations. 
In our investigation of (e + — H2) scattering, we considered the generalization of the 
Kohn method due to Kato [H EEO] as well as the complex Kohn method [TTJ [12] . Both of 
these modifications are designed to circumvent problems associated with Schwartz-type 
anomalous behaviour; the latter especially has become popular in recent years as it has 
been commonly believed automatically to be free of all nonphysical anomalies. 

We have carried out Kohn calculations for (e + — H2) scattering in the fixed- nuclei 
approximation [131 E] so that only the motion of the three light particles need be 
considered. All of our calculations use Hartree atomic units. Moreover, we regard the 
positron as being of sufficiently low energy that elastic scattering is the only significant 
open channel other than positron annihilation and that only the lowest partial wave 
of symmetry need be investigated. Armour and coworkers [15] have shown that 
consideration of this partial wave alone is sufficient to account for (e + — H 2 ) scattering 
processes up to ~ 2 eV. This partial wave is equivalent to the s-wave in atomic scattering. 

In this article, we will construct a framework ultimately to describe a formal 
connection between the mechanics of the real-generalized and complex Kohn methods. 
We will prove that the phase shift approximation determined in the complex Kohn 
method is exactly equal to that obtained from the real-generalized Kohn method 
(hereafter referred to as the generalized Kohn method, for brevity) under a particular 
optimization. In developing the framework leading to this conclusion, we will establish 
three other important results. Firstly, we will identify the mathematical origin of a 
certain class of anomalies as they appear in generalized Kohn calculations, before proving 
that such anomalies are necessarily absent in applications of the complex Kohn method. 
Secondly, we will show that the existence of the anomaly-free singularities we reported 
earlier [1] is fully accounted for by this framework. Finally, for both the generalized and 
complex Kohn methods, we will demonstrate that explicit solutions of the variational 
equations are not required to obtain optimized phase shift approximations, since these 
can be determined generally from the evaluation of four matrix determinants. The 
optimized results show that the generalized and complex Kohn methods are equivalent. 

The treatment given here for (e + — H 2 ) scattering is sufficiently general that it 
should not be difficult to adapt it to other physical systems. Indeed, it applies directly 
to potential scattering problems and single open channel scattering problems, with 
exchange if required. We will, therefore, concentrate on abstract mechanisms rather than 
presenting numerical data from individual calculations. Following the analysis, we will 
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briefly relate our framework to earlier work on avoiding the Schwartz singularities and 
the more familiar notation generally used by other authors (see, for example, Burke [T6] , 
and Burke and Joachain [17], with multichannel extensions summarized by Nesbet |18j). 
We will also demonstrate why we do not use this more familiar notation throughout. 

2. Anomalous behaviour in the generalized Kohn method 

2.1. Variational approximations to the phase shift 

Concerning the present calculations, our implementation of the generalized Kohn 
method is analogous to that described in a previous article PQ, so only a brief review 
of the basic principles will be given here. We approximate the exact leptonic scattering 
wavefunction, with a trial wavefunction, viz. 

M 

V t = (S + a t C)tp G + Y,PiXi, (1) 

i=l 

where 



s ' 




cos(r) sin(r) 




" S ' 


c 




— sin(r) cos(t) 




C 



(2) 



As we have explained elsewhere pQ, r G [0, n) is a phase parameter that can be adjusted 
to avoid nonphysical, anomalous behaviour in the results of Kohn calculations. The 
functions, S and C, form a basis that represents very low energy incident and scattered 
positrons asymptotically far from the target molecule; they account only for the lowest 
partial wave of E+ symmetry and they are independent of r. For convenience, their 
explicit forms will here be taken to be the same as those defined by equations (3) and 
(4) of our earlier article [I]. In those calculations, prolate spheroidal coordinates were 
used for their convenience in describing the two-centre H 2 system, though the general 
features of the following analysis are not particular to this choice. That is to say, without 
any great modifications an equivalent treatment could be given for different bases and 
coordinate systems better suited to other problems. For example, calculations on short- 
range radial-potential partial wave scattering [T6l [T7] use spherical coordinates with 
open channel radial functions which are asymptotically proportional to spherical Bessel 
functions. 

The function, ipQ, is a real- valued, unit-normalized approximation to the electronic 
ground state wavefunction of the unperturbed hydrogen molecule. It is evaluated 
at a fixed internuclear separation, typically the equilibrium value. Other than the 
requirements that it is square-integrable and independent of both r and the coordinates 
of the positron, the precise form of ipQ is not important here (although an accurate 
description of the target is, of course, important when calculating physically relevant 
results [19J). The functions, {xi\, are real-valued short-range correlation functions. 
They are used to describe interactions between the positron and the target electrons 
when they are close together. Apart from the restrictions that they must be independent 
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of r and square-integrable, thus contributing negligibly to fy t when the positron is 
asymptotically far from the target molecule, the precise forms of these functions are 
unimportant here. It should be noted that (Q]) does not explicitly contain the function, 
Xoj which appeared in the Kohn trial function used earlier pp. This omission does 
not affect the generality of our calculations and has been made purely to simplify the 
notation involved in the analysis that follows. 

As before [I], for brevity we will abbreviate by (X, Y ) any integral of the form 
(X\ (^H — E^j \ Y), where H is the nonrelativistic Hamiltonian for the scattering system 
in the adiabatic nuclei approximation and E is the sum of the kinetic energy of the 
positron and the ground state energy expectation value of tpc- Integrals of this form 
are evaluated over the configuration space of the positron and the two electrons. The 
operator, (jk — E^j, contains a term of the form — ~ (V 2 + k 2 ), where — |V 2 and k > 
are, respectively, the kinetic energy operator and the wavenumber of the positron. 
Consequently, for S and C as in pQ, it is well known that 

(S^ G ,C^ G ) - (C^ G ,S^ G ) = k, (3) 

where we have defined 

~ ttN 2 R 2 

k = ^f-k, (4) 

in which R is the fixed internuclear separation and N is an arbitrary normalization 
constant appearing as a factor in the explicit forms of S and C. The explicit form of 
the coefficient of proportionality in (j3J) is specific to the choice of basis. For example, 
in the aforementioned case of short-range potential scattering [TH1 [T7] the coefficient is 
simply N 2 /2, assuming the radial S and C functions are each associated with a common 
spherical harmonic for angle-function normalization. 

Since H is Hermitian, the {xi} are all square-integrable and each term in ([T|) is 
real-valued, we note that 



(Sip G ,Xi) =(Xi,Si/>G) (z = l,...,M), (5a) 
(CVG^HOfcC^G) (z = l,...,M), (56) 
(XhXj) ={Xj,Xi) (iJ = 1,...,M) . (5c) 

Henceforth, we will make implicit use of these properties. 

In the generalized Kohn method, a stationary principle is used to determine optimal 
values of the unknown parameters, a t and {p±, . . . ,Pm}, in ^t- This principle manifests 
itself in the system of linear equations, 

Ax = -6, (6) 



where 
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A 



x 



(Wg,W g ) (Ci/j g , X i) 
(Xi.CVg) (XiiXi) 

{Xm, Cip G ) (xm, Xi) 
(Ci/> Q ,Si(> G ) 

{XM,Slp G ) 

a t 
Pi 

Pm 



(Cipc, Xm) 
(Xi,Xm) 

\Xm, Xm) _ 



(7a) 



(76) 



(7c) 



If A is nonsingular then the solution of fl6]) uniquely determines optimal values for the 
unknown parameters in \l/ t . This solution can then be used to calculate a variational 
approximation, r] v G (— 7r/2, 7r/2], to the exact scattering phase shift, t]. We note 
for generality that if the method is applied to single channel electron scattering with 
exchange, whilst the definition of {X, Y) changes to account for antisymmetrization, the 
abstract forms of ()3]) and fl5aj) -f J7cj) are unaltered. 

For S and C as in [TJ, r] v is obtained implicitly from the definition, 



tan (r] v - r + c) = a t 



k 



(8) 



in which c = kR/2. The appearance of the parameter, c, in (|S]) is an artefact of the 
particular choice of basis functions, S and C, in our (e + — H2) calculations. When using 
prolate spheroidal coordinates, it is convenient to construct these functions in such a 
way that they include a phase factor equal to — c. The appearance of c in ([8]) merely 
compensates for this phase factor and plays no 'active' part in the following analysis. 
The functional, X [\l/ t ] , is defined as 

Z[tf t ] = (* t ,tt t ) (9) 

and, when the Kohn equations ([6]) are satisfied, can be shown to take the form 



M 



Z[* t ] = {8^ g Mg) +a t {Si> G ,Cip G ) +Y / Pj(^G,Xj)- 



(10) 



Under these circumstances, upon substitution of the solution of into ([8]), the error 
in tan (r] v — r + c) from tan (77 — r + c) can be shown to be second order in the error of 
\l/ t from \l/ [21 [91 [10] . This is the attraction of the Kohn variational principle. 
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2.2. Singularities and anomalous behaviour 

In the case where A is singular, the system of Kohn equations ([6|) either has no unique 
solution or no solution at all, and the variational method breaks down. Where A 
is nonsingular, it has been widely documented [TJ [Tl El EQl EU E21 [231 El] that the 
results of Kohn calculations can exhibit nonphysical behaviour when the values of 
various parameters in the Kohn trial function are close to values making A singular. 
In the generalized Kohn method, problems associated with these so-called Schwartz 
singularities are accounted for by the inclusion of the adjustable phase parameter, 
r G [0,7r). As we have already noted [I], for a sufficiently accurate trial function the 
choice of r should have no significant physical effect on the calculation of r) v and, in this 
respect, we can regard it as a free parameter. Consequently, if Ar] v denotes the difference 
in the phase shift approximations obtained at two similar values of r separated by At, 
the typical dependence of rj v on r is such that At] v /At is small. However, if r is close 
to a value, r s , making A singular, then nonphysical anomalies can arise and these are 
characterized by large values of Ar/ V / At. We have found P that a choice of r G [0, ir) 
can normally be made successfully to avoid anomalies of this kind at a given positron 
energy. 

We shall now develop a framework formally to describe the underlying mathematical 
structure of those anomalies in t] v (r) encountered by varying r at a fixed k. This 
framework offers a complete, analytic description of rj v (r). We begin by considering the 
determinant of A ( !7aj) ; from the Laplace expansion of det (A) along the first row of A , 
then subsequently along the first column of each resulting submatrix of A, using (j2J) we 
obtain 

M M 

det (A) = V (k) (C^g, C^g) + £ (*0 (Ctyo, Xi)(C^ G , Xj ) 

i=i j=i 

= A (k) sin 2 (r) + B (k) sin (r) cos (r) + C (k) cos 2 (r) , (11) 

where A (k), B (k), C(k),V (k) and Sy (k) = Sji (k) depend on k but are independent of 
t. If, by accident, k = k z such that A (k z ) = B (k z ) = C (k z ) = 0, then det (A; k = k z ) is 
identically zero and a unique value of r] v cannot be determined from the Kohn equations 
at any value of r. We define 

Z = {k z G M + : A (k z ) = B (k z ) = C (k z ) = 0}. (12) 

To ensure that our Kohn calculations can always be implemented consistently for 
some choice of r, we will henceforth consider only values of k G M + \ Z rather than 
k G M + , although for brevity we will not usually mention this distinction explicitly. In 
practice, we have not yet found any such value of k z in our numerical calculations. It is 
conceivable, though very unlikely, that such a value might occur by accident for some 
configuration of parameters in the trial wavefunction, but then it should be possible to 
vary the nonlinear parameters in the short-range correlation functions always to avoid 
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this eventuality. Hence, we consider the case k = k z here purely so that it can be 
excluded, in the interests of good mathematical practice. For each k G M + \ Z , there 
exists at least one value of r such that A is nonsingular. 

Next, consider the matrix, A, formed by replacing the first column of A by —b, 
where b is as in flTfel). so that 



A 



-(Xi,Sip G ) (xi,Xi) ••■ (Xi,Xm) 



(13) 



-{xm,Si[)g) (xm,Xi) ■■■ (xm,Xm) . 
By the same method used to find ffTTj) . we have 



M M 



det (ij = - V (k) (Ci/jq, Si/j g ) — ^2 Sij W (Ci>G, Xi) (SipG, Xj) 

i=i j=i 

= A (k) sin 2 (r) + B (k) sin (r) cos (r) + C (k) cos 2 (r) , (14) 

where A(k), B (k) and C (k) depend on k but are independent of r. The functions, 
V (k) and (k), are as in (fTTj) . Now, setting r = 0, r = 7r/2 and r = tt/4 successively 
in both (JTTi) and (114jl . it is readily shown that 



M M 

C(k) =V (k) (Cty> G , Cty G ) + Sij ( k ) X.) (C^g, Xi), 

M M 

i=l i=l 

M M 

^ (A;) = P (fc) (Sif> G , Si> G ) + W <^ G ' **> Wg, Xi), 

i=l i=l 

M M 

A(k) =V(k) (S^g,C^ g ) +J2J2 S ij ( k ) (S1>G,Xi)(C*l>G,Xi), 

i=l j=l 

B(k) = -P(fc)((5fc^ G ) + {^ G ,^ G )) 

M M 

~ ( k ) ({ClpG,Xi)(SlpG,Xj) + ( S ^G,Xi)(ClpG,X3)) 

i=l j=l 

B(k) =V(k)((SijG,S^ G )-(C^ G ,C^ G )) 

M M 

+ Y,ZZ S V ( k ) ((^G, Xi) Wo, X3) - (C^G, XiXC^G, Xi)) 

1=1 i=l 

from which we obtain the relations, 
A (jfc) - C (Aj) —B(k), 
A(k)-C(k) = -B(k), 



(15a) 
(156) 
(15c) 
(15(0 

(15e) 

(15/) 

(16a) 
(166) 
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so that we may rewrite (1141) as 

det (X) = A (k) sin 2 (r) + [A (k) - C (k)] sin (r) cos (r) 



+ 



A(k)+B (k) 



cos 2 (r) . (17) 
Now, for nonsingular A, we define (k) to be 

6 (k) = det (A) 1 [$? t ] . (18) 
Although det (A) and X [\&t] each depend on r, with effort it is possible to show that 



the product of these two terms does not. A proof of this result is given in Appendix A 
Next, by Cramer's rule, we can write 



det A 



det (A) ' 

for nonsingular A. Using f|T8|) and (j!9j) . for nonsingular A we can then rewrite (JS} as 

i -r(fc) 



(19) 



tan (rj v — r + c) 



det (A) 

where we have defined the function, T (k), which is independent of r, as 

Q(k) 



T(k) 



k 



(20) 



(21) 



Since, for each k G M + \ Z, it is always possible to find a nonsingular A for some choice 
of r, G (k) and T (k) are defined completely in this domain. Further, for S and C as in 
PQ, both G (k) and T (k) are continuous functions of k over R + \ Z. 
We will find it convenient to define the functions, 



f(r;k) = det(Aj -T(k) 

= A (k) sin 2 (r) + [A (k) - C (k)} sin (r) cos (r) 
+ \A (k) + B (k)] cos 2 (r) - T (jfc) 



and 



(22) 



g (r; fc) = det (A) = A (Jfc) sin 2 (r) + B (Jfc) sin (r) cos (r) + C (k) cos 2 (r) , (23) 

where we have explicitly denoted the parametric dependence of / (r; k) and g (r; k) on 
fc. For nonsingular A, we then have 

/ (r; fc) 



tan (r/ v - r + c) 

Next, we see that 

f (r; jfc) = [A (k) - C (jfc)] cos (2r) - B (jfc) sin (2r) 



(24) 



(25) 
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and 



g ' ( r; k) = [A (k) - C (k)} sin (2r) + B (k) cos (2r) , 



(26) 



where, in both cases, the prime indicates partial differentiation with respect to r. 
Moreover, differentiating (12^1) with respect to r, we find 



sec 2 (r] v - t + c) 



dr 



- 1 



P (r; k) 
9 2 (r; fc) 



^7v 
<9r 



f 



g(r;k) f (r; k) - f (t; k) g' (t; k) 
9 2 (r; fc) 



so that 



% / 2 (r; fc) + (r; fc) + ^ (r; fc) f (r; k) - f (r; fc) ^ (r; fc) 



dr 

It can be shown using 



5 yi ^^ if j 



/ 2 (r; fc)+s 2 (r;fc) 
and (T2SD that 



(27) 



(2* 



/ 2 (r; k) + g 2 (r; fc) + g (r; fc) /' (r; fc) - / (r; k) g' (r; fc) 
^ (k) - r (ife)l U (ife) - r (ife) + B (Jfe)l + ^ (fc) C (A;) , 



(29) 



which is independent of r. We define 

Q (k) = \A (k) - T (k)] \A (k) -T{k) + B (Jfe)l + ^ (fc) C (fc) (30) 



and, further, prove in Appendix B that 



G (k) = r 2 (k) . 

It follows immediately that 

dn v T 2 (k) 



(31) 



(32) 



dr Z 2 (t; *) + (r; *) ' 
for nonsingular A. 

Anomalous behaviour in t/ v (r) will arise whenever r is such that the ratio given by 
([32]) is unusually large. Singularities in a Kohn calculation appear whenever g (r; k) = 0, 
though it is clear from (l32il that the presence of singularities is neither a sufficient nor a 
necessary condition for anomalies to occur. For a fixed value of k, values of dr\ v jdr might 
be large even if g (r; k) is never zero. Conversely, for fixed k, large values of dr] v /dT are 
not guaranteed in the limit as r varies so that g (r; k) — > 0, owing to the fact that both 
/ (r; AO and g (r; k) appear in the denominator of the expression for dr) v /dr and their 
zeros will not coincide in general. This is the mathematical origin of the anomaly-free 
singularities reported previously PQ; in our own calculations, we have confirmed that 
those singularities are of such a type that, in the limit as r varies so that g (r; A;) — > 0, 
P (r; A:) does not become small in comparison to T 2 (k). 
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An immediate consequence of (1ST]) is that, using (I3"U|) and provided k is such that 
2A (k) + B (k) 7^ 0, we can write 



r(*) 



A(k) + B(k) A(k)+A(k)C(k) 



(33) 



2.A (k) + B (k) 

so that the functions of k in ( 1221) can be expressed purely in terms of ^4 (k), B (k), C (k) 
and A (k). Hence, in calculating the value of r) v from ( 124)) . in general there is actually no 
need to solve the Kohn equations fl6]). All that is required to determine r] v (r) completely 
at each k is to find the values of the four determinants, A (k), B(k), C (k) and A(k); 
substituting ( 133]) into ( 122|) . we find that ( 1241) can be rewritten 

[A (k) - C (k)} sin (r) cos (r) + B (k) cos 2 (r) + P (k) 



tan (?7 V 



r + c 



.A (k) sin 2 (r) + i3 (k) sin (r) cos (r) 
for nonsingular A, where we have defined 

A 2 (k) - A (k) C (k) 



C (k) cos 2 (r) 



If k = k s 
A(k 8 )- 



VVk) 

2A(k)+B(k) 

l M, + \Z such that 2 A ( k s ) + B {k s ) 
B(k s )]A(k s )+A(k s )C(k s ) =0. 



(34) 



(35) 



= 0, then from (IB.ll) and (IB.2I) we also have 
Here, the continuity of T (k) over R + \ Z 



ensures that the value of T (k) determined from (|33|) in the limit as k — > k s is well 
defined and equal to the value of T (k s ) determined using ( 12~TT) . ( 1A.8I) . ( 1A.17I) . ( 1A.18I) 
and (1A.20j) . In practical calculations, T (k s ) may be determined either by this method 
or by interpolation of T (k) either side of k s . The physical significance of instances of k s 
is discussed in section HI 



2. 3. Optimization 

Having found an analytical form for dr] v /dT, an obvious extension to our investigation 
is to optimize this expression with respect to r. From inspection of (|32|) we have 
dri v /dT > so, at each k, finding a global minimum of dr] v /dr with respect to r G [0, 7r) 
locates the point at which r] v (r) varies most slowly with r. This forms a natural 
optimization scheme for choosing r to avoid anomalous behaviour. 

For nonsingular A, partial differentiation of ( 1321) with respect to r gives 

ffin Y 2 (k\ 

Or 2 [p ( T ; k) + g 2 (r; k)] 

and, after some manipulation, we find 

/ (r; k) f (r; k) + g (r; k) g' (r; k) = X (k) sin (2r) + y (k) cos (2r) , (37) 
where we have defined 

x[k) = AHk)-B>(k)-cHk )+B{k) 



r (k) - A(k) 



(38) 
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and 

y (k) = [r (k) - A (Jfe)l [C (k) - A (k)} + A{k)B (k) , (39) 

both of which are independent of r. Denoting by Tj any value of r for which 
d 2 rj v /dT 2 = 0, if r (k) ^0 we then have 

X (k) sin (2n) + y (k) cos (27*) = 0. (40) 

In the special case where k = k g such that r (fc g ) = 0, then dr) v /dr is everywhere 
zero and optimization is not required since rj v is constant over r. In the special case 
where k = k&, such that X (fc h ) = y (k h ) = 0, then for nonsingular A the value of 
d 2 rj v /dT 2 is everywhere zero. Hence, dr\ v /dr is constant with respect to variations in r. 
If this constant is equal to zero, optimization of dt] v /dT is not required since then t] v is 
constant over r. If this constant is nonzero, rj v varies linearly with r and no preferred 
optimization can reasonably be defined. 

Discounting these two special cases, at each k there will be exactly two values of 
Tj G [0, 7r) satisfying ( I4"0"l) . separated by n/2. In this case, using ( 1371) and ( HOI and 
differentiating (1361) with respect to r, at r = Tj and for nonsingular A we see that 

_ , dr2m [* (fc) cos (2rQ - y (fc) sin (2rQ] 
^ 3 (T " Tt) " W [/Mr^)+. 2 (r^)] 2 ' ( } 

In general, when k is such that both T (k) and at least one of X (k) or y (k) is nonzero, 
then (HOI) ensures that d 3 r) v /dr 3 is nonzero at r = Tj. Moreover, since the two values of 
T are separated by n/2, we see from (HI!) that the signs of d 3 rj v /dT 3 at the two values 
of Tj are opposite. Hence, in general, drj v /dT has one minimum and one maximum 
for t G [0, 7r). At each k, we will denote by To and T\ the values of Tj respectively 
minimizing and maximizing d7] v /dr. We will denote the values of n v (T = r ) and 
rj v (t = Ti) respectively by 774° and r/i 1 ' 1 . 

Next, assuming that A is nonsingular and k is such that To and Ti exist, we note 
the following. Firstly, we see from ( 138]) and ( )39|) that these conditions preclude having 
k = k c such that both ^4 (A; c ) = C (fc c ) and B (k c ) = 0, since X (k c ) = y (k c ) = by 
inspection. Consequently, we note from ( |25|) and ( |26i) that /' (t; fc) and g' (t; fc) cannot 
both be zero. In fact, f'(rf,k) cannot be zero, since this would require g'(Tf,k) ^ 
and, using ( )40l) . we see that ( 1371) could not then be satisfied, since we have assumed 
that g (n; k) ^ 0. Hence, using (j2i|) . (1371) and (]30]1 . we write 

/' (T i5 k) tan (r/« - t + c) + g 1 (t; fc) = 0, (42) 

where /' (t; fc) ^ 0. Both t and r± must satisfy ( 1421) . Since To and T\ are separated by 
7r/2, using (I25p and (T26l) we can then immediately conclude 

tan (r/( 0) - t + c) = tan - n + c) , (43) 

so that, for r/ v G (— t/2, t/2], the values of and t]^ 1 '' must also be separated by n/2. 
In fact, using (1251), (1251), ij3"gjh (}39lh (TO and (1321), it is straightforward to show that 
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w _ , , [A{k)-C{k)]y{k)-B{k)X{k) 



tan ( V ; Tl + c) [A (jfe) _ c ( ^ + ^ ( fc) y (jfc) 
2,4 (*;)+£ (A;) -2r(fc) 



*4(fc) + C(fc) 
3. The complex Kohn method 

5. i. Variational approximations to the phase shift 

The complex Kohn method is an extension of the original variational approach in 
which the boundary conditions of the trial wavefunction are complex-valued. Although 
originally thought [ITJ [12] entirely to be free of Schwartz-type behaviour, anomalies 
have been reported by Lucchese [25] and, recently, by Cooper and coworkers [1]. We 
now consider an implementation of the complex Kohn method which is analogous to 
that described previously pQ. Following the same approach as in the preceding section, 
we now develop a mathematical formalism to explain the success of the complex Kohn 
method in avoiding a particular class of anomalies. The same notes and minor caveats 
concerning applications of the method to systems other than (e + — H2) apply as in 
subsection 12.11 

We begin with a complex- valued trial wavefunction, viz. 

M 

* t = {S + a[T)^ G + J2p^ (45) 



i=i 



where 



T = S + iC, (46) 

the functions, S, C, ipQ and {xi, ■ ■ ■ , Xm} being as in (JT]). The unknowns, 
{a[,p[, . . . ,p'm}j wm n °t) m general, be real. The primes on {a' t ,p[, . . . ,p' M } distinguish 
them from the corresponding values found in ([T]). Throughout this section, unless 
otherwise noted we will use primes in this way to distinguish various quantities used in 
our application of the complex Kohn method from the corresponding quantities used in 
our application of the generalized Kohn method involving a real- valued trial function. 

Application of the Kohn stationary principle leads to a set of equations analogous 
to ©, 

A'x' = -b', (47) 

where 



A' 



(f^ G ,f^ G ) (T*Vg,Xi) 
(Xi,Tip G ) (xi,Xi) 



(XM,T?p G ) (xm,Xi) 



(T^G,XM> 

(Xi,Xm) 



(48 a) 



(Xm, Xm) 
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b' 



(xm, Si> G ) 



(486) 



x 



Pi 



I Pm 



(48 c) 



Here, T* is the complex conjugate of T. In the usual Dirac notation, (T| implies complex 
conjugation of T. However, as pointed out by Chamberlain [26] (see also [HI [12]), in a 
consistent implementation of the complex Kohn method this conjugation of the 'radial' 
function should not, in fact, be performed. Hence, we have replaced (T\ by (T*\ to 
indicate that the conjugation is not carried out. 

If A' is nonsingular then the solution of (|4T|) uniquely determines optimal values for 
the unknown parameters in This solution can then be used to calculate a variational 
approximation, 77^ G C, to the exact scattering phase shift. For S and C as in [I], this 
estimate is obtained implicitly from the definition, 



tan (r]' v - r + c) 





% 




1 + a[ + ik- 1 !' 





(49) 



where 



X' 



(** t \(H-E)\* t ) = {*l* t ) 



(50) 



is analogous to (Q. In subsection 13.21 we demonstrate that the denominator of (1491) 
is nonzero when both A and A' are nonsingular. When the Kohn equations ( 1471) are 



satisfied, X' 



X' 



takes a form analogous to ( flOl) . 



M 



(51) 



Under these circumstances, upon substitution of the solution of fj47j) into fT49l) . the error 
in tan (77^ — r + c) from tan (77 — r + c) can be shown to be second order in the error of 
^ t from 

It should be noted that the value of r]' Y will not, in general, be real. However, since r\ 
must be real, the imaginary part of rj' v can be regarded as an error term arising from the 
fact that the trial function, is inexact. We will identify its precise form in subsection 

El 
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3.2. Avoidance of anomalous behaviour 

In the case where A' is singular, the system of Kohn equations ( )47|) either has no 
unique solution or no solution at all, and the variational method breaks down. For 
nonsingular A', we now demonstrate that the value of r]' v obtained in our implementation 
of the complex Kohn method is independent of the choice of r. We consider first the 
determinant of A'. Proceeding in a manner analogous to subsection I2.2[ it is easily 
shown that 

det (A') = [A (k) -C(k)- iB (k)] exp (-2ir) , (52) 

so that, as noted before [I], det (A') describes a circle in the complex plane for variations 
of r G [0, 7r). Hence, in the complex Kohn method, singularities can neither be located 
nor avoided by varying only r. Here, A (k), B (k) and C (k) are as in ( fill . 

Next, consider the matrix, A', formed by replacing the first column of A' by —b', 
so that 



A' 



-(T^ g ,^g) CT*V>g,Xi) ••• (T*V>g,Xm) 
-(XuStpc) (Xi.Xi) •■■ (Xi,Xm) 

-(Xm,Si/>g) (Xm,Xi) ■■■ (xm,Xm) 



(53) 



Following the same approach taken in subsection I2.2[ after a little work it is possible to 
show that 



det [ A' 



iA (k) -C{k) + [C (k) - A (k) + iB (k)] cos (r) exp (-ir) 



(54) 



where A(k), B (k) and C (k) are as in (TIT]) and A (k) is as in (T141) . Next, we will find it 
convenient to define the functions 

u (r; k) = -iC (k) - A(k) + [iC (k) - iA (k) - B (k)] cos (r) exp (-ir) + T (k) (55) 
and 

v (r; k) = [A (k) - C (k) - iB (k)] exp (-2ir) - iu (r; k) , (56) 
noting that these functions satisfy the identity 

u 2 (r; k) + v 2 (r; k) + v (r; k) v! (r; k) — u (r; k) v' (r; k) = 0, (57) 
where the primes on u' (r; k) and v' (r; k) indicate partial differentiation with respect to 



r. Next, in Appendix A for nonsingular A' we prove that 



det (A')l' 



(5« 



where 9 (k) is as defined in (fT8|) . Using Cramer's rule, together with (l2Tj) . (152|) . (154|) . 
fl55l). fl56l) and fl58l). we find that fT49l) can be rewritten 



tan (?7y - r + c) 



idet i' +T(fc) 



u (r; A;) 



det (A') + det (i') - ir (k) v ( T > k ) ' 



(59) 
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provided that v (r; k) is nonzero. We note in passing that inspection of (13"3"j) . (15 2p . (154"]) 
and fl59|) indicates that the phase shift approximation in the complex Kohn method 
can be evaluated generally from the same four determinants required in the case of the 
generalized Kohn method and without the need to solve the Kohn equations (|47|) . 
By analogy with ( 1281) . we then see that 



dy' v = u 2 (r; k) + v 2 (r; k) + v (r; k) u' (r; k) - u (t; k) v' (r; k) 

Or ~ u 2 (r;k) + v 2 (r;k) ' [ ' 

Inspection of fl52l) and fl56l) shows that the zeros of u (r; /c) and v (r; A;) coincide if and 
only if A' is singular. Hence, for nonsingular A', the denominator of flBT?]) is nonzero so 
that, using ( 157|) . (160]) becomes 

^ - 0, (61) 



giving 



<9r 



^ ^ °- (62) 

Thus, whenever det (A') and v (r; /;;) are both nonzero, the value of i]^ is independent 
of the choice of r in ( 1451) . Complex Kohn calculations of 77^ will automatically 
be free of those Schwartz-type anomalies characterized by large values of dr\ N jdr in 
the implementation of the generalized Kohn method already discussed in section [2j 
Nevertheless, anomalies could still arise in the results of complex Kohn calculations due 
to the choice of some other parameter in the trial function; it is likely that this is the 
underlying cause of the persistent anomalies described earlier pQ. 

If v (r; k) = 0, then from fl57j) we either have u (r; k) = v (r; k) = 0, in which case 
det (A') = 0, or we have u (r; k) = v' (r; k). We also find 

3 [v 1 (r; k) - u (r; k)\ = det (A) , (63) 

where A is as in fTTajl . Hence, v (r; k) is zero only if at least one of A or A' is singular. We 
will shortly carry out a formal comparison of the results of the generalized and complex 
Kohn methods. In so doing, we will take the parameters of our trial functions to be 
such that both A and A' are nonsingular, so that the Kohn equations in each case can 
uniquely be solved. These conditions automatically ensure that v (r; k) ^ 0. 



3.3. Equivalence 

We now demonstrate the effective equivalence of the generalized and complex Kohn 
variational methods. Consider the two trial wavefunctions, \l/ t and which contain 
the same approximate target wavefunction and identical sets of short-range correlation 
functions. Suppose that A (JToj) and A 1 (148 ah are nonsingular so that solutions of the 
Kohn equations, (Q and f|4"Tj) . uniquely exist and we have v (r; k) ^ 0. Suppose further 
that k ^ k g and k 7^ k^ in order that ( l32i) can uniquely be minimized, as discussed 
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in subsection 12.31 These conditions are sufficient for (j44p and ( 15 9 p to be well defined. 
Under these circumstances, we claim that 

^W v ]=vi° ] - (64) 

Here, we can regard 9ft [r) ; v ] as the approximation to the phase shift in the complex Kohn 
method, since we have already noted that the imaginary part of rj^ can be interpreted as 
an error term. We recall that rj^ is that value of r] v obtained in the generalized Kohn 
method at the unique value, r = To, which minimizes dr] v /dr. 

Proof. We consider tan (rj® -rf Y -n + A. Using (|22), 422), (SID, 455), (ESJ) and (EH}, 
together with the standard result 

tan(jP Q) ~ l + tan(P)tan(Q)' (65) 
after some considerable manipulation, it can be shown that 

»[tan(#-<-r 4 + r)]=^ (66) 



and 



(67) 



where we have defined 

a (r; fc) = # (fc) sin (2r) + y (Jfe) cos (2r) (68) 

and 

6(r;A;) = -/ 2 (r;A;)-/(r;A ; ), (69) 

noting that b (r; fc) < 0, since we have assumed g (r; fc) 7^ so that A is nonsingular. 
Derivations of (166]) and (IBTl) are outlined in Appendix C Now, setting r = Tj, we have 
a{ji]k) =0 from ( HOI . Hence, 

sin - 3? [7^]) cos (r/^ - 3? [/?;]) 
» [tan (rtf> - V ' v )] = V ? \ V L = (70) 



cos 2 (r$> - 3? + sinh 2 (9 [ V ' v ]) 



and 



3 [tan (r/W - V >)] = sinh (9 ^]) cosh(9 [<]) = < Q 

cos 2 (r# } - [r/;]) + sinh 2 (3 [t£]) 6 ( r ^ fc ) 



(71) 

Taking G (-7r/2,7r/2] and 3* [<] G (-vr/2, vr/2], since 77^ and rj^ are separated 
by 7r/2 we can immediately conclude from (17D1 that we have either 3? [77^] = t]^ '' or 
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3ft [rj' v ] = rji 1 ^. Moreover, using ( l32l) and (l69~l) . by the definitions of To and t\ it is plain 
that 

T 2 (k) T 2 (k) . . 

< , , —k < °> ( 72 ) 



b(rr,k) b(r ;k) 

noting from (1701) and (ITTj) that 6 (r ; k) and 6 (ri; fc) cannot be equal, since rff 1 and rji 1 ^ 
are separated by n/2 and exactly one of r) v °^ or 77^ must give cos (^qi^ — 3ft [r/^] J = 0. 

It follows directly from (17T1) and (1721) that cos 2 ^77^ — 3ft [rj' y ]\ > 0. Consequently, 
inspection of (1701) reveals 



3ft [r}' v ] = n 



(o) 

V 5 



as required. □ 

An interesting consequence of this result is that, setting r = 7$ and then i = in 
(IPTj) . using 032]), ([64]) and ([69]), it is evident that 

tanh(SK]) = ^(r = r ), (73) 



so that we may write 

~dr) v 



ril 0) +itanh _1 



dr 



to) 



(74) 



Hence, the imaginary part of the complex- valued approximation to the scattering phase 
shift obtained in the complex Kohn method can be used as a measure of the susceptibility 
of the corresponding generalized Kohn calculation to Schwartz-type behaviour. In the 
case of a calculation involving the exact scattering wavefunction, the imaginary part of 
r]' v would be zero and the corresponding generalized Kohn calculation of rj v would be 
independent of r. 



4. The relationship to previous work on the generalized Kohn method 

In this section, we briefly relate our current work to earlier studies of Kohn-type methods 
and Schwartz singularities [U [TJ El EDI EU [221 E31 El]. This is intended as a guide and we 
do not attempt the full rigorous approach of the rest of the article. For clarity, we use 
notation similar to that of Burke and Joachain [TTJ. The multichannel extension of this 
notation is given, for example, by Nesbet [18] and Lucchese [25]. Note that in [17] the 
operator, 2 yE — H^j , for a short-range radial potential scattering problem is considered 

rather than (^H — E^j as used here. For completeness, we note that alternative versions 
of the Kohn method have been developed in terms of a Feshbach projection operator 
formalism [27] and have been found [28] to give anomaly-free results. More recently these 
methods have been revived and further developed [29] and have produced accurate phase 
shifts for low energy electron hydrogen atom scattering [30]. 
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The traditional approach is to separate the 'closed' part of the matrix, A, from 
the 'open' part which comprises the first row and column and involves the functions, 



S and C. This may be done by inverting the 'closed' matrix (Arl in Appendix A ) 



or, equivalently, by diagonalizing it. This separation aids the analysis but may not 
necessarily be carried out in practical calculations as each term in A is energy- dependent, 
although the inverse or diagonalization need only be calculated once for all r. The closed 
terms are then 'folded' into the open channel matrix elements as optical potentials 
[TBI IT?] [18] so that the generalized Kohn problem is reduced to one involving matrices 



of a dimension equal to the number of open channels. If the eigenvalues of A^l are 
taken as (Ej — E), f = 1, . . . , M, then we define 

t ici q i \ ST (^g,X/)(X/,^g) 

L n = (SiffQ, SipG) - 2_, ^ _ e) ' ( ' 

^22 = (Ci/} G , Cty> G ) - 2^ ^ _ ^ , (756) 

^12 = (Stye, CVg) - 2^ _ ^ , (75c) 

^21 = (C^G, Stye) - 2^ (jjT _ ) > ( 75fi? ) 

where the are the diagonalized linear combinations of the %%■ 

In an obvious notation, the are related to their r = values by 



(76) 





L\2 


H 


\ L 21 


L22 J 






and the generalized Kohn result is 



tan (r/ v — r + c) = — -=^i — ^ — det (L) , (77) 
£22 kL 22 



in which 

det (L) = L n L 22 - L 12 L 21 (78) 

is independent of r from ( 1761) . We note that (L i2 — Z21) and (Z 2 2 + -^11) are also 
independent of r. 

The problem with this approach is that it introduces singularities at the energies, 
E = Ef, which need to be accounted for. In practical calculations, E will certainly 
range across one or more of these poles as the variational open channel functions are 
required to take into account only the asymptotic behaviour of the exact open channel 
functions (although more sophisticated functions than this may be chosen). When the 



Equivalence of the generalized and complex Kohn variational methods 



19 



Kohn method was first introduced these poles were briefly considered as causes of the 
Schwartz-type anomalies. However, Nesbet [HI [20] showed that (I77|) is nonsingular at 
these energies; the second order poles in det (L) cancel leaving first order poles which 
cancel with those in the denominator. 

The framework of this article avoids the universal introduction of these poles and 
also avoids the need for closed channel diagonalization to relate the analysis to results 
of practical calculations. Intermediate poles are limited to expressions that use (133]) to 
determine T (k s ), and we have described a nonsingular expression for T (k s ) at the end 
of subsection I2.2[ using the work of Appendix A 



We may illustrate the relationship between the current approach and the traditional 
approach by considering an expansion of det (A) with the closed matrix, A^l, replaced 
by the diagonal matrix of eigenvalues, viz. 



M 



det (A) =(C4j G ,C4i G )l[(E f -E) 



M M 
/=1 /'=!,/¥/ 



or, introducing the intermediate poles, 

~ M 

n ( E f - e ) 

L/=i 



det (A) 



L 



22- 



Forming a similar expression for the matrix, A, and writing 



M 

F (*) = II (E f -E) = det 



/=i 



we have the correspondence, 



(79) 



(80) 



(81) 



Also, 



A (k) ^F(k)L u , 

C(k) ^F(k)L 22 , 

B(k) <— > -F(k)(L 12 + L 21 ) 

A(k) ^F(k)L 12 . 



2A (k) + B{k)^F (k) (L 12 - L 21 ) = F (k) k, 

so that each instance of k = k s corresponds to E = Ef for some /, and 

det (L) 



T(k) <—>F(k) 



k 



12 a) 

m) 

32 c) 

12 d) 

(83) 
(84) 
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We have reproduced all of the results for the generalized and complex Kohn methods 
derived in the preceding sections independently using the formalism, though not 
always with the same strict rigour. We note that the relatively laborious proofs in 
Appendices A and B are required for strict avoidance of intermediate poles, otherwise 
the algebra is of equivalent complexity. We note that in the case of the more familiar T- 
matrix and .S-matrix versions of the complex Kohn method, with open channel functions 

S + aJ(C + iS) (85) 

and 

(C-iS) -af(C + i§), (86) 

respectively, the behaviour with respect to r is the same as presented here, with correct 
to second order variational estimates 

al = 1 + 2iaJ (87) 

and 

aZ = i(a'Y 



as expected. Here, aj and af play a role analogous to a' t in (Hoi) , and we have defined 



a' v = a[ + il' 
k 



(89) 



This leads to the T- matrix method variational estimate of the phase shift, 77J, having 
the form 

r /v T = r / (°)-itanh- 1 ^(r = r )l. (90) 

We end this section with a few remarks on the earlier ways the generalized Kohn 
method was used to avoid anomalous behaviour [201 1221 [231 121], as summarized by 
Nesbet [18] . These various methods certainly avoid the anomalous singularities but 
they generally do so by attempting to maximize the absolute value of det (A) or L22 as 
a function of r. We may take the derivative of (1TTT) with respect to r and set it to zero. 
Denoting by r d any value making det (A) stationary with respect to r, we then obtain 
the expression 

[A (k) - C (k)] sin (2r d ) + B (k) cos (2r d ) = 0. (91) 

Thus, if k ^ k^, there are exactly two values of r G [0, 7r) making det (A) stationary, 
separated by n/2. We label these values and Td 2 , and note that they are distinct 
from t and T\. We abbreviate by and rj^ 2 ^ the values of 77 v (T dl ) and rj v (r d2 ) 

respectively. In general, including the particular case where the two stationary values 
of det (A) are of equal magnitude and opposite sign, 77V and r/i^ are not equal (see 
[3T] . which contains material supplementary to that presented here and in [I]). Certain 
of the anomaly-free methods have developed additional techniques to cope with this 
possibility pi 
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Remembering that det (A) is periodic in r, the conditions for stationary values 
with respect to r are related to the conditions for existence of the Kohn singularities. 
The two stationary values of det (A) will, in general, comprise one maximum and one 
minimum. These extrema will either have the same sign (no singularities), or exactly one 
of them will be zero (exactly one singularity), or they will have different signs (exactly 
two singularities). This is obvious qualitatively, and a detailed algebraic description 
of the behaviour of det (A) with respect to r is available from the authors. From the 
argument we have proposed [1] regarding the existence of anomaly-free singularities, 
in the case above where exactly one singularity exists, for a sufficiently accurate trial 
wavefunction we would expect it to correspond to an anomaly-free calculation of the 
phase shift approximation. In the case above where exactly two singularities exist, for 
a sufficiently accurate trial wavefunction we would expect one of these singularities to 
be anomaly- free and the other to be anomalous, in the sense discussed earlier pQ. 

Some of the above referenced anomaly-free methods attempt to use r as a 
variational parameter once the anomalous region has been avoided. This is done either to 
help increase the phase shift, though we note that the Kohn methods provide stationary 
but not bounded variational estimates, or otherwise to try to improve the 'quality' of the 
trial wavefunction according to various criteria put forward. This is in contrast to the 
current work, in which we argue that r should ideally be an arbitrary parameter and we 
look for the minimum value of dr] v /dr. The 'quality' of the scattering wavefunction then 
depends on the forms of the S, C and {xi} functions when the projectile is close to the 
target. We suggest that, in addition to the direct relationship with the complex Kohn 
method derived above, improvements in computing power and computational science 
over the last few decades justify this approach, with nonlinear parameters in the bulk 
trial functions available for use as additional variational parameters. 

Cooper [31] has carried out some studies of the behaviour of dr] v /dr at r = r , as 
a function of k. This behaviour is related to the separation of the anomaly-free and 
anomalous singularities as a function of k [H E]. As defined [31], this separation is a 
value in the range [0,7r/2]. When the two singularities are separated by 7r/2, their r 
values satisfy the condition (HP]) required of r and T\. However, dr] v /dr at r = r is 
generally smaller as a function of k when the separation of the two singularities is also 
small [3T] . 

We note that the earlier anomaly-free methods have been extended to the 
multichannel case [T8 l 1231 1251 132] . The complex Kohn method is being applied to various 
calculations of electron polyatomic molecule scattering (see, for example, [331 [M]). We 
are currently extending the framework and analysis of this article to the multichannel 
case and we aim to develop the analysis of the anomaly-free singularities as part of the 
multichannel work. 
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5. Concluding remarks 

We have shown that, in the case of phase shift calculations for low energy (e + — H2) 
scattering, the complex Kohn method is equivalent to a particular optimization of the 
generalized Kohn method. Further, we have established a number of interesting results 
regarding the appearance of anomalous behaviour in our generalized Kohn calculations. 
Specifically, we have found that anomalies which appear when only r is varied can be 
explained from purely analytic considerations; they are intrinsic to the Kohn method 
itself and do not, as has previously been suggested [35J, arise from matrix ill-conditioning 
or errors due to limited computational precision. Our analysis describes analytically the 
behaviour of the phase shift over the entire range of r. This makes it possible to give 
a full description, for the first time, of any anomalous behaviour that results from the 
variation of r. 

By obtaining an analytic expression for drj v /dT, we have explained the 
mathematical origin of the anomaly- free singularities identified in our earlier article [I]. 
This result complements the physical argument for the existence of these singularities 
given in that article. 

We have demonstrated that there is a particular class of anomalies that are 
necessarily avoided in our complex Kohn calculations. However, it is important to 
note that this method, as we have implemented it here, avoids only those anomalies 
encountered in generalized Kohn calculations by varying r and keeping other free 
parameters, such as k, R and the nonlinear parameters in the trial function, fixed. 
Anomalous results characterized by unusually large values of d^R.[r]' v ] /dk, say, could 
still arise even in the complex Kohn method. We have not developed an explicit 
expression for this derivative in the way that we have here for d7] v /dr and dr]' v /dr; 
strictly speaking, we have not even considered in detail whether <99ft [rj' v ] jdk exists in 
the sense of demonstrating that 5i [rj v ] is a differentiable function of k. However, were an 
analytic expression for d^t [rf y ] jdk available, it is reasonable to conclude that it would 
explain the persistent anomalous behaviour discussed in pQ. 
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Appendix A. The function 

For nonsingular A, we claim 



det {A)X[%]=Q{k^) t 




where we have used a self-evident notation to denote that (k,X) is independent of r. 
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Proof. In the following argument, we will write = O (k, r) before explicitly proving 
the independence of from r. Henceforth, it will be convenient to consider only the 
case M > 3 in (CQ). It is straightforward to show that the result ( {TBI) is satisfied for 
M < 3 by explicitly inverting the matrix, A, allowing X[\l/ t ] to be found. Throughout, 
we will implicitly make use of the Hermiticity properties ( l5aj) -(!5cj). 

We denote by the (M + 1 x M + 1) matrix formed by replacing the j th column 
of A f Tffflj) by —b ( 176)) . We will denote by Aj?l the (M x M) matrix formed by removing 
the i th row and j th column of A. The row and column indices of A range from 1 to 
M + 1. By assumption, A is nonsingular, so that using ( ITU]) and Cramer's rule, which 
states that 

det f A 



a ' = (A - 2) 



and, for 1 < j < M, 



det ( 



» = (A - 3) 

the product, det (A)Z[\l/ t ], can be written 



(k, r) = det (A) 1 [* t ] = det (A) (Sty Gj ^g) + det (A (1) J (S^ G , Cty G ) 

+ det (^O+i)) (^G,Xi>- (A.4) 
i=i 

The Laplace expansion of det (A) along column 1 of A is 

M 

det (A) = (Cfa, C^ G ) det + ]T (C^, x.) det (a^) , (A.5) 

i=l 

while the expansion of det yA(j)j along column j of Ayj is 

M 

(Ctyo, det (^gj) + ]T (-1/ (^ G , X.) det (a[ 
Now, using (T5j), we obtain 



det ( A(j) 



i=i 



. (A.6) 



(Si(> G ,Sil> G )(CiJ> G ,Cil> G ) - (Sip G ,Ci; G )(C^ G ,S^ G ) 
= (Sij G ,Sij G )(Ctfj Gl Ctlj G } - (^ g ,CVg)(C^g,^g), (A.7) 

which is independent of r. Hence, by noting that det (^(i)J is a l so independent of r, 
then defining 
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x det [A 



and combining ( 1A.4I) . ( 1A.5I) . ( 1A.6I) and ( 1A.7I) . we can write 



M 

6 (A;, r) - 6o (k,S) = {Sip G , S^ G ) £ (-1) 4 (CV>G, Xi> det (a£ 



(i+l) 
) 



i=l 



- (SipQ, Cip G ) ("I/ (5Vg, Xi> det ( A 



1=1 



J'=l 



(^ G ,^G>det Ug^ 



+ ^(-ir(^ G , X ,)det(A (i+i n 



cj+i) y 



(A.9) 



It then remains to be shown that the right hand side of (1A.9I) is independent of r. 

Using the fact that det (M) = det (M T ) for any square matrix, M, together with 

the fact that A is symmetric, we deduce that det f A^j^ = det • We can then 



rewrite flA.91) as 



' M 

E 

i=i 



-lVdet U 



I (i+l) 
L (l) 



x ((Sil> G ,Sif> G )(Cif> G ,Xi) 
- (S^ G ,Cip G )(Sip G ,Xi) 

(Ci> G ,Si> G )(Si; G ,Xi)) 



M M 

, ^i+j+l ( (i+l) 



i=l j=l 



(A.10) 



We next define the (M x M) matrix, X = A^l. Further, we denote by xfy the 
(M — 1 x M — 1) matrix formed by removing the z th row and j' th column of X. Further, 
for i ^ p and j ^ q, we denote by X^ p } the (M — 2 x M — 2) matrix formed by removing 
the i th and p th rows and j th and q th columns of X. The row and column indices of X 
range from 1 to M. The elements of X, xQ and X?j p } are independent of r. If % = p 

or j = q, we define X^'^ to be the (M — 2 x M — 2) matrix of zeros. We find that 



detU (m) ^ ' ^ W 



l (i) 



i=i 



(A.n; 
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det (Ag$) = det (xgj) (CV G ,Cy G > 



M M 
_p=l g=l 



P+Q+1 + 



x detfX^J(^ , % ,><C r Vo,X,> 



where terms of the form have the definition 



(A.12) 



0~ab 




1 



(a > 6) 
(a < b) 



Using (lA.lip and f !A.12j) . we rewrite flA.lOp as 

e (fc, r) = e (k,r) + e x r) + e 2 r) 

where 



(A.13) 



(A.14) 



M M / 

9 1 (A;,r) = ^^(-ir +J+1 det(x ( ( i]) x (% ) ^ G )(^ G)Xi )(^ G)Xi ) 

t=i j=i V 

- (Si/j G ,Ci/j G ){Stp G ,Xi)(CiljG,Xj) 

- (C , ^G,5 , ^ G )(5'^G,Xi)(C'^G,Xj) 



and 



e 2 (fc,r) 



M M M M 
i=l jr' = l p=l q=l 



0.9) 



X (^G, Xi){SlpG, Xj)(C*p G , X P )(ClpG, Xq) 



We will now show that 0i (/c, r) and 9 2 (fc, r) are each independent of r. 



Considering (IA.15j) . since X is symmetric, we note that det (X^M = det (x$ 



(A.16) 



Using (J2J), when the summation over % and j in flA.15j) is carried out, a number of terms 
cancel. We then find 
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M M 



e 1 (^ = ^^(-i) i+J+1 det(x«) x Us^s^ G )(c^ Gl Xi)(c^Xj) 



i=i j=i 



(Si/; G , Cip G )(Sij G , Xi) (CiI>g, Xj) 

(Ci/jq, Sifj G ) (Sip G , Xi) {Cip G , Xj) 
(Cip G ,Cip G ){S^ G ,Xi)(S^ G ,Xj) 1 , (A.17) 



which is independent of r. Applying the same method to (1A. 16j) . with a little work we 
can write 



M M M M 
i=l j = l p=l q=l 



Cm) 



X (^G, Xi) (S^G, Xj) (CV'G, Xp) (Clp G , Xq) 



(A.18) 



which is independent of r. The cancellation in the summation ( 1A.16I) arises from the 
fact that 



+ (-l) CT ^ =0 (a ^ b) . 
Finally, combining ( 1A.14|) . ( 1A.17I) and ( 1A.18j) . we have 

e (k,y) = e (k,*) + e x (k,S) + e 2 , 

so that = 6 (k,x), as required. 



(A.19) 

(A.20) 
□ 



In the case of the complex Kohn method, a result similar to ffTBl can be derived by 
a method analogous to that given above. For nonsingular A', we claim 



det (A')X' 



-e . 



Proof. We define a function, A (k, r), such that 



A(fc,r) =det (A')l' 



(A.21) 



(A.22) 



Next, after some manipulation, it is straightforward to show that 



(^g,^g)(T> g ,T^g) - (^G,T^ G )(f> G ,^ G ) 
= (^g,^g)Wg,^g) - (^ g ,^g)<^g,C^g) 
= (Sty G ,Cty G )<Cty G ,Sty> G ) - (Sfc, Sip G )(C^ G , CiJ> G ), 



(A.23) 
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where we have used (IA.7j) . Proceeding in a manner analogous to that used above and 
adopting the same notation, using flA.8j) it is clear that we can then immediately write 



where 



A (k, t) = -9 (k,r) + Ai (k, t) + A 2 (k, r) 



(A.24) 



Ai(fc,r) 



M M 



^^(-l) iW1 det(x|j x ( ^ G) %>(TV G , Xi )(r> GlXj > 
i=i j=i 

- (T*i(> G ,Sil> G )(Sfc,Xi)(T*ii>G,X3) 



and 



A 2 (fc,r) 



M M M M 

j=l j = l p=l (?=1 



x (SipG, Xi) (Sipc, Xj) (r*Va, Xp) ( t >g, Xg) 



(A.26) 



Ax (k, t) and A 2 (k, r) being analogous to ( IA.15j) and (1A. 16j) . respectively. Considering 
first Ai (k, r), we have 



(S2pG,SlpG)(T*iJG,Xi)(T*^G,Xj) - (S^G,Ttp G )(SiJ G ,Xi)(T*^G,Xj) 

- (f*fc,S1> G )(Si(> G ,Xi)(f*il>G,X3) + (f*il> G ,Tii> G )(Sil>G,xd(Si!>G,Xj) 

= (S^ G ,Clp G }(SlpG,Xi)(C^G,Xj) - (S^G,Sip G )(Cip G ,Xi}(C^G,Xj) 

- (Cjp G ,CjpG)(SiJG,Xi)(s^G,Xj) + (c^s^ G )(s^x i )(CiJG,x 1 ) 

+ i{Sij> G ,Si!> G ) ({C^G,Xi){S^G,Xj) - (S^ G ,Xi){C^G,Xj)) ■ (A.27) 

When the summation in (1A.25I) is carried out, it is clear that the final terms in the 
square brackets in ( 1A.27I) sum to zero. Using (1A.15I) and (1A.17I) . we then have 

AI AI I 

A 1 (^ = -^^(-l) l+J+1 det(x ( ( ; ) ) ) x (S^G, S^ G )(C^ G , Xi)(C^G, Xj) 

i=i j=i V 

- (SlpG, Clp G ) (SlpG, Xi) (ClpG, Xj) 

- {C^G,Sifj G }(Sil} G7 Xi){CipG,Xj) 

+ (CiPG,Cip G )(S^G,Xi)(SipG,Xj)) , (A.28) 
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A 1 (k,S) = -Qi (k,S) . 

Next, we consider A2 (k,r). Clearly, 



2n 



(A.29) 



(Sip G , Xi) (SipG, Xj) (T*ip G , X P ) (T*ip G , Xq) 
= (SipG, Xi) (SipG, Xj) (SipG, X P ) (SipG, Xq) 

- (SipG, Xi) (S4>G, Xj) (ClpG, Xp) (ClpG, Xq) 
+ i(StpG, Xi) {SipG, Xj) (SipG, Xp) (ClpG, Xq) 

+ i(SipG, Xi) (SipG, Xj) (Cip G , Xp) (SipG, Xq) ■ 



(A.30) 



Using ( 1A.19I) . we see that the first term in the expansion ( 1A.30[) does not give an overall 
contribution to the sum (1A.26I) . For the same reason, each of the final two terms in 
(1A.30I) also sums to zero in ( 1A.26I) . Hence, using ( 1A.16I) and ( 1A.18I) . we then have 



M M M M 
i=l j = l p=l q=l 



X (SipG, Xi) (Sip G , Xj) (ClpG, Xp) (Cip G , Xq) 



so that 



and, finally, 



as required. 



a 2 (k,S) = -e 2 (k,s> 



A(k,s) = det(A')!' % 



-6 (k,?) 



(A.31) 

(A.32) 

(A.33) 
□ 



Appendix B. The Q function 

For T (k) and Q (k) as defined in (1211) and f)30p respectively, we claim that (13T|) holds. 
Proof. To prove (|3T|) . inspection of ( )30|) shows that it is sufficient to prove 

H = 0, (B.l) 
where we have defined 



n 



A{k) + B{k) A{k)+ A{k)C{k)- 2A(k) + B(k) T (k) . (B.2) 



Henceforth, it will be convenient to consider only the case M > 3 in (TjQ). It can be 
shown that the result fIB.ll) is satisfied for M < 3 by explicitly evaluating expressions 
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for A(k), A (k), B (k), C (k) and T (k). Throughout, we will implicitly make use of the 
Hermiticity properties fl5aj) - fl5cj) . 

First, using f ll5al) . f ll5cl) . ( 115dl) and (11 5 el) , we obtain 



A (k) C (k) = V 2 (k) (SipQ, S^ G ) (Ctf> G , C^ G ) 

M M 

+ V (k) (S^g, S^ G ) J2J2 S v ^ (We* Xi)(CipG, Xj) 
i=l j=l 

M M 

+ V (k) (CVg, Cty G ) E E S V ( fc ) <^ G > Xi> (^g, Xj) 

i=l j=l 

M M M M 

+ EEEE^ ( k ) S pq( k ) (S*pG,Xt}(SlpG,Xj}{ClpG,Xp}{C^G,Xq} 
i=l j = l p=l q=l 



(B-3) 



and 



M M 



A(k) + B (k) = -V (k) (Ctyo, StPg) "EE S V ( fc ) (^g, X*> (StPg, Xj)- 

i=t j=l 



Next, in the nomenclature of Appendix A we note that 
V (k) = det (A$\ = det (X) 

and 



S ij {k)=S ji (k) = (-l) i+i+1 det[X i 



-0 

-U) 



1 



1 det(x$ 



(B.4) 



(B.5) 



(B.6) 



Combining (llSdjl . flB~3|) and (JR4]), as well as using f lA~8]) . flA~T7j) . flB~5|) and ( IB~6|) . after 
a little work we can write 



A(k) + B (k)\ A (k) + A (k) C{k)=V (k) [9 (fc) + 6i (Jfe)] (B.7) 

M M M M 

+ E E E E 5 « (*o 5 w (*o ^^ g ' **) ^^ g ' ^^ g ' x ^ ^^ g ' *«) 

8=1 j = l p=l g =l 
M M M M 

~ E E E E Si i W S P1 ^ ^ G ' ( C ^G, Xj) (C^G, Xp) (SlpGi Xq)- 
i=l j = l p=l q=l 



Next, summing (\15d\i and (IB. 41) . we find 
2A (k) + B(k) — V (k) k, 



(B.8) 



where we have used (j3J), together with the fact that Sij (k) = Sji (k). Recalling ( 12TT) . it 
is clear that 



2A(k) + B(k) r (k) = V (k) 9 (k) 



(B.9) 
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Using (IA.18p . (IA.20p . (IB.7P and (IB.9p . after some cancellation we can then rewrite (IB. 21) 

as 

M M M M 

H = EEEE T (-l) i+j+P+q {Srfa, x.) (Srfa, Xj ) (CVg, Xv) (Cifa, Xq ), (B. 10) 

i=l j'=l p=l g=l 

where, using ( IB .5 1) and ( IB. 6ft . we have defined 



T = det ( Xg) det - det (X$ ) det 



HO 



- (-1)^+^ det (X) det (Xgjj) • (B.ll) 



If i = p or j = q, it follows trivially from the definition of -Xy'gj given in Appendix A 



that T = 0. When i < p and j < q, we make use of the following result, 

det (X) det (Jtgg) = det (xgj) det - det det (*g>) . (B.12) 

Commonly known as the Lewis Carroll identity after its role in Dodgson condensation 
[36] , Fomin and Zelevinsky [37] point out that (IB. 12|) was, in fact, proved earlier by 



Desnanot (see, for example, [38]). It is easily seen to generalize to any % ^ p and j ^ q 
by multiplying the left hand side of (1B.12P by a factor of (— 1)°> +Cr ^_ Hence, from 
inspection of ( IB. lip we see that T is identically zero and, since 

T=0=>H = 0=>G(k) =T 2 (fc), (B.13) 

the required result (}3~Tj) follows. □ 

Appendix C. Derivations of (1661) and (167ft 

We consider a calculation under the same conditions on k as outlined at the beginning 
of subsection 13.31 Under these circumstances, we claim that ( 16B1 and ( 1671) hold. 
To derive ( ]66|) . we begin by defining the following functions, 

V (k) =2A{k) + B(k)-2T(k), (C.l) 

£(k) =A{k)+C{k), (C.2) 

p (r; k) =V{k)v (r; k) - E (k) u (r; k) , (C.3) 

q ( T ; k) = V (k) u (r; k) + £ (k) v (r; k) , (C.4) 

where u(r;k) and v(r;k) are as in (l55|) and ( 1561 . respectively. Using OH]), (JSHJ) and 
(I65p . we then find 

R [tan f„«> - ^ - r + r)l - — — ^ — — - — — + ^ ^ fel] _ ^ fell f C 5) 

L K ^ v ri + rjJ -R[g(r;fc)]R[g(r;fc)] + S[g(r;fc)]S[g(r;A;)]- 

We further define 
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d 1 (r; k) = [A (k) - C (k)] cos (2r) - B (k) sin (2r) 
d 2 (r; jfc) = [A (k) - C (k)] sin (2r) + B (jfc) cos (2r) 

After a little work, we then find 



(C.6) 
(C.7) 



3* b (r; fc)] K [g (r; fc)] + 3 [p (r; fc)] 3 [g (r; fc)] 
[D 2 (ifc) - S 2 (Jfc)] [3ft [u (r; jfc)] di (r; jfc) - 9 [u (r; jfc)] d 2 (r; jfc)] 

£> (Jfc) £ (fc) [d 2 (r; jfc) + d 2 (r; jfc) + 23 [u (r; jfc)] d x (r; jfc) + 23ft [« (r; jfc)] d 2 (r; jfc)] (C. 



and 



3ft [g (r; fc)] 3ft [g (t; fc)] + 3 [q (r; fc)] 3 [g (r; fc)] 
= 2V (k) S (k) [3ft [u (r; jfc)] d x (r; k) - 3 [u (r; jfc)] d 2 (r; jfc)] 
+ S 2 (k) [d\ (r; jfc) + d\ (r; jfc) + 23 [u (r; jfc)] di (r; jfc) + 23? [u (r; jfc)] d 2 (r; jfc)] 
+ [V 2 (k) + £ 2 (k)] [3ft [« (r; fc)] 3ft [« (r; fc)] + 3 [u (r; fc)] 3 [u (r; fc)]] , 

where we have made use of (1561) to eliminate t> (r; k). Next, we define 



(C.9) 



M (k) = A (k) - T (jfc) [C (k) - A (k)] +B{k)C (k) , 

M(k) =B(k)\A(k)-r(k)] + A2 W + B2 W- c2 W 



(CIO) 
(C.ll) 



so that, after some laborious but elementary operations, using f|55|) we obtain 

3ft [u (r; jfc)] di (r; jfc) - 3 [u (r; fe)] d 2 (r; k) = M (k) cos (2r) + A/" (ife) sin (2r) , (C.12) 

together with 



d 2 (r; jfc) + d 2 (r; jfc) + 23 [u (r; jfc)] d x (r; jfc) + 23ft [u (r; jfc)] d 2 (r; jfc) 
2A4 (jfc) sin (2r) - 2JV (fc) cos (2r) 



(C.13) 



and 



3ft [u (r; jfc)] 3ft [u (r; jfc)] + 3 [u (r; k)\ 3 [u (r; jfc)] = 2Af (k) cos 2 (r) - M (k) sin (2r) 



+ 



A(k)-T(k) +C 2 {k) 



Next, we define 



Z (jfc) = -4r 2 (jfc) - B 2 (jfc) - [„4 (jfc) - C (A;)] 2 , 
so that, using {33]), (EHJ), ([39]) and ffOTSl) - ( f07T5|) . after a lengthy effort we find 



(C.14) 



(C.15) 
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M [p (r; k)} & [g (r; fc)] + 3 [p (r; fc)] 3 [g (r; fc)] 
Z (fc) [AT (k) sin (2r) + y (k) cos (2r)] 



and 



(C.16) 



» [g (r; fc)] 3ft [g (r; fc)] + 9 [g (r; fc)] 9f [g (r; 



2# (jfe) cos 2 (r) - y (fc) sin (2r) - i (ifc) - T (k) - A 2 (k) 



(C.17) 



Finally, after some work, examination of ( 122]) and (123]) yields, 



-f (r; jfe) -# 2 (r; jfe) = 2* (ifc) cos 2 (r) - J> (jfe) sin (2r) - ^ (ifc) - T (jfe) -^ 2 (jfe) .(C.18 



The required result (J55) then follows immediately from (155}. (155]). flCl5]) . (1U17]) . (ITHSl) 
and cancellation of Z (k). Since we have assumed k ^ k g and k ^ fch, it follows directly 
from flSE]), (ES} and flCTol) that £ (jfe) < 0. 

The result (167]) is derived by first noting that 

3 [tan f„W - r/ - r + r^l - ^ Ip (r; k)}^ (r; k)} - [q (r; k)Wp (r; k)} 

[ [Vv Vv Tl + rjJ " K [g (r; ifc)] K [g (r; jfe)] + 9 [g (r; ifc)] 9 [g (r; Jfc)] • (U1Jj 

By a method analogous to that used to derive (1ST)]) , we eventually find that 

9 [p (r; fc)] & [g (r; fc)] - 3 [g (r; fc)] & [p (r; fc)] = Z (A;) T 2 (k) . (C.20) 
The required result ( I6T}) then follows in the obvious way. 
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